Response in kinetic Ising model to oscillating magnetic fields 
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Ising models obeying Glauber dynamics in a temporally oscillating magnetic field are analyzed. In 
the context of stochastic resonance, the response in the magnetization is calculated by means of both 
a mean-field theory with linear-response approximation, and the time-dependent Ginzburg-Landau 
equation. Analytic results for the temperature and frequency dependent response, including the 
resonance temperature, compare favorably with simulation data. 
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I. INTRODUCTION 

Ising models with Glauber dynamics in an oscillating 
magnetic field were recently considered with Monte Carlo 
(MC) simulations in jfij and (2) . The phenomenon of 
stochastic resonance (pTwas revealed, exhibiting a char- 
acteristic peak in the correlation function C(T) between 
the external oscillating magnetic field and the magne- 
tization M(t) versus the temperature T of the system. 
The resonance temperature T r (the temperature at which 
C(T) has a maximum) was systematically computed as 
a function of the driving period, lattice size and driving 
amplitude, both for two-dimensional (2D) |J and three- 
dimensional (3D) |Q systems. The one-dimensional (ID) 
case was analysed by Brey and Prados Q in a linear-field 
approximation. 

The present work is a natural continuation of those 
studies, considering analytically the 2D and 3D cases. 
We will present two approaches. The mean- field the- 
ory with linear response approximation will be discussed 
first. Then in 2D where the mean-field theory is not 
as good as in other dimensions, a more refined time- 
dependent Ginzburg-Landau (TDGL) approach will be 
presented, with significant improvements. 

Recently, kinetic Ising systems in oscillating external 
fields have also been examined both experimentally and 
theoretically in 0. The focus was on properties below 
the zero-field critical point, such as the frequency depen- 
dence of the probability distributions for the hysteresis- 
loop area and the residence time. The latter quantity for 
small systems in moderately weak fields suggests further 
evidences of stochastic resonance. Very recently, finite- 
size effects versus driving frequency are analyzed as a 
dynamical critical phenomena ||. In contrast to these 
works, ours is focused on the temperature dependence 
above the zero-field critical point. 



II. MEAN-FIELD THEORY AND 
LINEAR-RESPONSE APPROXIMATION 

Our starting point is the master equation for the ki- 
netic Ising model obeying Glauber dynamics : 

P(a;t+1)- P(a;t) = 

]T [w(a' -> a)P(a'; t) - w(a - a')P(a; t)] , (1) 

a' 

where P(<t; t) is the joint probability of finding the spin 
configuration a at time t, and w's are the transition rates 
between two configurations which differ by one spin flip. 
For the heat-bath algorithm, the rate function is chosen 
as 



w(o — ► a ) 



1 



1 + e -/3[E(<r)-E{<7>)] ' 



with j3 = 1/T (hereafter the Boltzmann constant k = 1), 
and E(a) is the energy of a in a magnetic field h: 



E(a) = -jJ2SiSj-h(t)J2Si 



(2) 



where h(t) — Asin(u)t) and X)nn denotes a summation 
over nearest neighbors in a square or cubic lattice. 

Let us denote the configuration a by the values of the 
spins Si, Si, Sy, with system volume given by V = 
N d . d is the spatial dimension of the system and TV is its 
linear size. Since Si = ±1, it is easy to rewrite (^) as 

, v 
j t P{S u S 2 , S v ;t) = -J2^j(Sj)P(Si,S2,...,S v ;t) 



3=1 



^2 W 3(- S j) P ( S l, S 2, ■■; ~Sj, S V ; t) 



(3) 



with 



^j(S j ) = ^[l-S^ a xih(E j /T)}, 



Ej = jJ2S k + h, 



(4) 



fe=i 



1 



where the last sum runs over the z nearest neighbors of 
the spin Sj, with z — 2d. Multiplying both sides of (||) 
by Si and performing an ensemble average (denoted by 
< • • • >), after some simple mathematical tricks, we get 
the basic equation for the Glauber dynamics: 



dt 



< S, >-- 



< Si> + < tanh(Ei/T) > 



(5) 



Invoking the mean-field approximation, we replace Ei by 
Jz < S > +h to get: 



dt 



<S>= - <S> +tanh[(/i + < S >)/T], (6) 



where T C MF = Jz is the mean-field critical temperature. 
In the absence of h, the magnetization is given by the 
stationary solution of the well-known equation: 



< S > = tanh[T c MF < S > /T] 



(7) 



For small h(t), we may use the linear-response theory 
in (§) by first writing < S > (t) =< S > +AS(t) 
and considering the h/T -C 1 and AS/T -C 1 limits. 
Performing the Taylor expansion and keeping only the 
first-order terms, equation (0) becomes 



dt 



AS 



+ -(l-<S> 2 )sm(u J t), 

T~MF J 



where 



tmf 



(l-<5>g) 



(8) 



(9) 



is the relaxation time. The solution can be found easily: 



AS(t) = AS sin(u>t - MF ), 
with the phase shift and amplitude given by 
#mf = arctan(wTMF) 
AS = ^(1- < S > 2 ) 



1 

MF 



(10) 

(11) 
(12) 



The correlation function between the total magnetiza- 
tion M = V < S > and the external field h(t) can be 
computed: 



C = M(t)h(t) = (Vlu /2tt) 
VA 2 



2T 



■(!-<£>$) 



2tt/u 


tmf 



AS(t)h(t)dt 



1 + wV, 2 



(13) 



MF 



Here the overline denotes a temporal average over a pe- 
riod P = 2tt/lj. In the T > T C MF domain, < S > = 0, 
thus C becomes: 



VA 2 



T -T, 



MF 



2 (T - T C MF ) 2 



2^2 ' 



(14) 



III. TIME-DEPENDENT GINZBURG-LANDAU 
APPROACH 



Before comparing (|13|) to simulations, we present an 
alternative, continuum approach to compute C. For 
an Ising system with non-conservative order parame- 
ter (model A Q), the time-dependent Ginzburg-Landau 
(TDGL) equation for the local magnetization density 
<p(f, t) takes the following form: 



90 

at 



9 jA 



(15) 
(16) 



where ri is the coarse grained Hamiltonian. For our 
present purpose, the white noise £(r, t) which accounts 
for the effect of thermal fluctuations is irrelevant. Con- 
ventionally, parameters T, u and g in ( |l6| ) are understood 
to be obtained by coarse graining the microscopic dynam- 
ics ([!]). For critical properties, the important tempera- 
ture dependence in these parameters lies inuoc T — T C GL , 
giving rise to the spontaneous symmetry breaking be- 
low the critical temperature T C GL . In order to compare 
with simulations, more precise dependences on T are re- 
quired. To this end, we outline here a refined mean-field 
approach in the continuum limit. The same approach 
has been successfully applied to the two-species driven 
diffusive systems This approximation is expected to 
be good outside the critical region. However, this turns 
out to be not a serious handicap because the presence of 
an oscillating field prevents the system from building up 
critical correlations. 

In a mean-field approximation, the joint probabilities 
in (|l|) are factorized into singlet probabilities p(r; t) for 
finding the spin up at site r at time t. This effectively 
produces the power series expansion of ri in 4>. 1 — p gives 
the probability of finding the spin down. The continuum 
limit involves an expansion in the derivatives, such as: 



p(x±l,y;t) ^p(x,y;t)± 



dp(x,y;t) 1 d 2 p(x,y;t) 



Ox 



dx 2 



By identifying p as (0+ 1) /2, we obtain from (|l|) a kinetic 
equation for 4> after some algebra. For h = 0, we find 
precisely the deterministic part of ([l5]) with: 

r = h-2W i + 2W- i -W 8 + W- 8 ), (17) 
8 

u = -^(6VF + 12W 4 - 4VF_4 + 5H^ 8 - 3VF_ 8 ), (18) 
81 



2r 



(-6Wo-4W 4 + 4W_4 + 5W8 + PF-8), (19) 



where W n = 1/(1 + e nl3J ) contains the desired explicit 
T dependence. When a small uniform field h is applied, 
to 0(h) (neglecting a <f> 5 term) we have finally the deter- 
ministic kinetic equation 



2 



dt X 



6 



(20) 



where /x = /3(3VF 2 + 4W 4 W_ 4 + W 8 W_ 8 )/2T. It is useful 
to note that T, 5 and \x in (|2^) are positive definite for all 
T, whereas u has one zero at T C GL « 3.0901J w 1.3618T C , 
where T c = -2/ln(\/2 - I) J w 2.2692 J is exact. This 
is an improvement over F = 4 J from the last section. 
Moreover, we reproduce the first few terms of the high- 
temperature series expansions of thermodynamic quan- 
tities such as the susceptibility and the relaxation time. 
In the f3 — ► limit, we recover the mean-field results of 
the last section: u « 1//3J — 4, T w ft J, g w 48(/?J) 2 , 
and /x w 1/ J. 

For small h and T > T C GL , the nonlinear term (j</> 3 
in ( po|) is negligible. The total magnetization M(t) = 
/ drtfi(r, i) — <p(q — 0, t) in response to an external field 
can then be computed easily, where 4> denotes the spatial 
Fourier transform of <fi. It satisfies dM/dt = — TuM + 
Tfih(q = Q,t). We readily find 



M(i) 



VfiAT 



■.sin(ut- 6» G l), 



(21) 



where the phase shift is 6*gl = arctan(w/r-it). The corre- 
lation function with h is then given by 



Ct>t gl — 



VA 2 T 2 iiu 
2[{Yu) 2 + oj 2 



(22) 



Note that this coincides with the mean-field result (|14| ) 
in the high-temperature limit. 

For T < T C GL , the term proportional to g is needed 
to break the sym metry, le ading to the spontaneous mag- 
netization m = \J —6u/g. Linearizing about m, we find 
precisely the same form of C as T > T C GL except that u 
is replaced by — 2u in (22). 



IV. DISCUSSION AND COMPARISON WITH 
SIMULATIONS 

^From the simulation data in |jj and j^j , we learn that 
the system has a maximum response to external driving 
at a definite temperature T r which depends on the driving 
frequency. Hence T r can be designated as the resonance 
temperature. ^From the analytically determined correla- 
tion functions in ( |l4| ) and (|||) , we find two peaks in C at 
above and below the respective T c , and also C(T C ) = 0, 
as shown in Fig. 1. This double-peak structure in C is 
consistent with simulations for larger lattice sizes (up to 
N = 200 for 2D and N = 40 for 3D) and with smaller 
steps in T than reported in jij and The reason for 
missing the peak below T c in our earlier simulations is 
due to the use of small lattice sizes. Note that the peak 
below T c is much smaller than the one above and its po- 
sition is less sensitive to the driving period. The reason 
for the overestimated theoretical values of the peaks be- 
low T c may be accounted for by the frustration of the 



system to order in the presence of hit); such frustration 
steming from metastability has not been taken into ac- 
count explicitly in our theories below T c when we seek 
the symmetry-breaking solutions. 

We believe that this also explains the discrepancy at 
T c , where simulations show a small but finite C(T). 
Finite-size effects are not of great concern here because, 
as mentioned above, the correlation length even at T c 
is truncated by h. In simulations, we have checked the 
convergence in C(T) for N > 50 in 2D. 

Focusing on T > T c from now on, the TDGL pre- 
dictions for C(T) are more accurate than those of the 
mean-field theory in general. They both converge to the 
simulations in the tails at T > T c (see Fig. 1). In 3D the 
mean-field theory is already acceptable. 

Turning our attention to the amplitude dependence, 
replotting the simulation data from JlJ and [gj suggests 
that the height of the peak C(T r ) cx A 2 , in agree- 



ment with Ql4) and (|22j). For not too large frequen- 
cies and small A, the theoretical proportionality constant 
agrees well with simulations. For example, the slope of 
C(T r )/(VT c ) versus A 2 /T 2 for P = 50 in 2D gives 0.92 
from simulations |5J], 0.96 from TDGL and 0.99 from 
mean-field approach. In 3D the same slope is 0.88 from 
simulations £3], and 1.29 from mean-field approach (In 
3D the comparison are worse because T r is much closer 
now to T c .) This proportionality is a manifestation of the 
linear response of the system to h, which breaks down at 
large enough amplitudes. Our new simulations show that 
this happens for A/T c > 0.15 in 2D for P = 40. 

A quantity of significant interest is the resonance tem- 
perature T r (P). It can be determined analytically from 



T 



MF 



T. 



MF 



1 



(23) 



and numerically from (p2|) for T GL . These together with 
simulation results are presented in Fig. 2. The agree- 
ments are reasonable. As expected the mean-field ap- 
proximation is quite good in 3D but in 2D the TDGL 
approximation is better. 

The results in Fig. 2 confirm the earlier observation in 
and U that for P -> 00 we get T r -^T C . This result is 
also consistent with the one obtained by Brey and Prados 
Q in ID where the above limit becomes T r — > T c = 0. In 
the opposite limit P — ► 1 (in unit of Monte Carlo steps 
P > 1) both the theory in ID and our approximations 
in 2D and 3D suggest T r — > const. Unfortunately, in jlj 
and Q the wrong conclusion T r — > 00 was drawn in this 
limit. Similarly, the position of the peak below T c also 
converges to T c in the P — > 00 limit. 

In passing, we also derive fiofl the relationship between 
the correlation function and the hysteresis-loop area A: 



A = 2itC I tanf 



(24) 



where 9 is the phase shift between h and M. This then 
relates our results of C to that of A as observed in || . 
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V. CONCLUSIONS 



Using mean-field with linear-response and TDGL ap- 
proximations, the characteristics of the resonance peaks 
observed in kinetic Ising models in oscillating magnetic 
fields are reproduced. New simulations improve ear- 
lier results by confirming the analytically predicted dou- 
ble peaks. Focusing mostly on the behavior above T c 
(where our approaches work better), we determine the 
dependence of the resonance temperature as a function 
of driving frequency and amplitude. We confirm the al- 
ready predicted result in |]|| that T r — > T c for the limit 
of practically interesting driving frequencies (P — > oo), 
and corrected the wrong extrapolation in the opposite 
limit P — > 1. We introduce a refined TDGL approach 
which improves significantly the mean- field results in 2D, 
but in 3D the mean-field approximation is already ac- 
ceptable. We have thus demonstrated that the stochastic 
resonance in kinetic Ising models above T c can be under- 
stood by means of rather simple theoretical approaches 
for small driving amplitudes. 
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FIG. 1. C(T)/(VT C ) versus temperature for P = 40 and A = 0.05T C for 2D in (a) and 3D in (b). Dots are MC simulation 
results in 2D (N = 200), triangles are MC simulations in 3D (N = 40), continuous line is from TDGL approximation and the 
dashed line is the mean-field result. The higher peak for TDGL than mean-field theory below T c is due to our dropping the (f) 5 
term in (20). 
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